Reconstructing cosmic acceleration from f(R) modified gravity 
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A variant of the accelerating cosmology reconstruction program is developed for f(R) gravity. 
Reconstruction schemes in terms of e-foldings and by using an auxiliary scalar field are developed. 
An example of a model with a transient phantom behavior without real matter is explicitly discussed 
in both schemes. Detailed comparison of the two schemes of reconstruction is presented. It seems 
to support physical non-equivalence of the two frames. 
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I. INTRODUCTION 
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One of the most important problems of modern cosmology is the explanation of the current universe speed-up, first 
discovered in [l|, Q . A convenient way to express this situation is to introduce a new form of energy, called Dark 
Energy (DE). During the last few years several theories have been developed in order to formulate and try to explain 
■ the dark energy universe. 

| As gravitational alternatives for DE, modified gravity theories have been formulated, calling for plausible late-time 
modification of General Relativity (GR). Many modified gravity models have been proposed in the literature (for a 
review, see Q), starting from the very simple 1/R theory 4-6] to more elaborated ones. 

For any such theory to be valid it is always strictly required that it accurately describes the known sequence of 
cosmological epochs, specifically it must accurately fit an increasing number of more and more precise observational 
data @,[8j . Another important issue to be taken into account concerns the different types of finite-time, future 
cosmological singularities (see @). It has been demonstrated (see [l(| El) that in modified gravities all four known 
types of finite-time future singularities may possibly appear. However, modified gravity can have the chance to cure 
qH' all of these future singularities (see [uH Il3j). 

Several reconstruction methods that make possible to reproduce any given cosmology have appeared for f(R) 
(2JT). modified gravity (l4T - [l^ | . In the following these methods will be shown. 

II. COSMOLOGICAL RECONSTRUCTION OF REALISTIC MODIFIED f(R) GRAVITIES IN TERMS 
^ ■ OF E-FOLDINGS 

In this section, the techniques of |l4j, [15| are summarized in order to construct an f(R) model realizing any given 
CNI ■ cosmology. The starting action for f(R) gravity (see e.g. Q, for a general review) is 

g i S = J d A x^—g + C matter \ . (1) 
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Using the e- folding variable, N = In assuming the matter density p is given in terms of a sum of fluid densities 
with constant EoS parameters, u>i, and considering the Hubble rate is given by G(N) = H 2 , the first Friedmann 
equation for the action ([T]) yields 

= -9G(N(R)) [4G'(N(R)) + G"{N(R))} + [3G(N(R)) + §G'(A(i?))] « - ^ + 

+ Eftoao 3(1+Ws) e- 3(1+Wl)7V(i?) - (2) 

i 

This is a differential equation for f(R), where the variable is the scalar curvature R, which is here R = 3G" (N)+12G(N) 
(for detailed calculations see 

Let us consider an evolution given by the following Hubble parameter: 

H 2 (N) = H Q In (^) + Hx = H a N + H X = G(N), (3) 
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where Hq and Hi are positive constants. We can, in this case, achieve a transient phantom behavior without the 
presence of real matter. Indeed, from R = 3G'(N) + 12G(N), we can obtain N — N(R) and Eq. (HH takes the form 
(for more details, see [l5[) 

= -3H (R 3ff )^ + (*=P + f) d Jm _ £M , (4 ) 
whose solutions are given by the Kummer's series $(a, b; z), namely 

f(R) = C<S (-2, -I; ^a) = C (-| + ^i? - ^i? 2 ) , (5) 

where C is a constant. As a consequence, with this f(R) theory, as given by Eq. ([5]), we can reproduce the transient 
phantom behavior without real matter given by Eq. © . 



III. COSMOLOGICAL RECONSTRUCTION OF REALISTIC MODIFIED f(R) GRAVITIES USING A 

SCALAR FIELD 



In this section it will be used the technique of [HI, [l6[ in order to construct an f(R) model realizing any given 
cosmology. We start from the action for modified gravity 



S ' = J d 4 Xy/^g(f(R) + Cmatter), 



which is equivalent to 



S = / d*Xs/^{P{4>)R + Q(^) + C matter ). 



(6) 



(7) 



Here, Cmatter is the matter Lagrangian density and P and Q are proper functions of the auxiliary scalar field, <f>. After 
some calculations and some assumptions (see fl6|) we can obtain the following differential equation for the function 

p(4>) 



-3(1+0)*) 3(1+« ( ) S (<A) 



(8) 



While the function Q{4>) is given by 



Q(4>) = -6(</(0)) 2 P(0) - 6. 9 '(0) 



3(l+aO 3(1+^)3(0) 



(9) 



The function f(R) can be calculated as f(R) = P(<j>(R))R + Q(<f)(R)). Thus, any given cosmology, expressed as 
a(t) = aoe 9 ^, can be realized by some specific /(i?)-gravity. 



We consider another time the transient phantom behavior, without real matter, given by ([3]). In this case, Eq.® 
reduces to 







d 2 P(4>) H , 



(10) 



We can obtain P(4>) from this equation and Q((f>) from ©.The function f(R) — P(4>(R))R + Q(<fi(R)) reduces to (for 
more details, see fla]): 



f(R) = Ci A + H oy /3{R-3y/2H ) 



El 
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A - El 

A 24 



3(H-3y'2H ) 



-R-3 % /2ff 



24-2(i?-3\/2Ho) 



2 V 6 V 2 



R-3,/2ff 



2 U 6 V 2 



z^R-SVZHo J e-y'dy), 



(11) 



/ 



where ^^S^y-^ + ^-J. We thus prove that, with this scheme, we are able to obtain the f(R) model (JTTTJ) which 
reproduces the desired transient phantom behavior without real matter, as given by ((31). 
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IV. SUMMARY AND DISCUSSION 



In summary, two different schemes of reconstruction have been developed for modified gravity. With these methods, 
any explicitly given cosmology can be realized as a corresponding modified gravity. 

The same example has been developed in the two cases. The result obtained in the first scheme (in terms of 
e- foldings) is given by Eq.© while for the second scheme (using an auxiliary scalar field), the f(R) obtained is given 
by (IIIII) . As one can easily see, the results obtained for both methods are in fact different. The reason behind this is 
the fact that action fljjl corresponds to a wider class of theories than the action ([6]) (for a related and quite detailed 
discussion, see [13, Ell). Nevertheless, if in Eq. we set C2 = 0, then the results coming from both schemes are 
similar, at least in the sense that, for low curvatures, they behave as constant, while for large curvatures the behavior 
is in both cases as R 2 . 

Finally, I want to remark that this work is based on the work done with Emilio Elizalde, which was published in 
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